Coherent phonon scattering effects on thermal transport in thin semiconductor 
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The thermal conductance by phonons of a quasi-one-dimensional solid with isotope or defect 
scattering is studied using the Landauer formalism for thermal transport. The conductance shows 
a crossover from localized to Ohmic behavior, just as for electrons, but the nature of this crossover 
is modified by derealization of phonons at low frequency. A scalable numerical transfer-matrix 
technique is developed and applied to model quasi-one-dimensional systems in order to confirm 
simple analytic predictions. We argue that existing thermal conductivity data on semiconductor 
nanowires, showing an unexpected linear dependence, can be understood through a model that 
combines incoherent surface scattering for short-wavelength phonons with nearly ballistic long- 
wavelength phonons. It is also found that even when strong phonon localization effects would be 
observed if defects are distributed throughout the wire, localization effects are much weaker when 
defects are localized at the boundary, as in current experiments. 

PACS numbers: 74.50+r 74.72-h 



I 

C 

O 
o 



> ■ 
■ 

o : 
r» . 

o ■ 

: 

o . 



i 

C 

O 

o 



X 



I. INTRODUCTION 

The transport of heat by phonons in quasi-one- 
dimensional structures is sensitive to quantum confine- 
ment effects once the phonon wavelengths that domi- 
nate thermal transport are comparable to the structure 
dimensions. Strong quantum confinement in thermal 
transport has been demonstrated in the observation^ of 
the thermal conductance quantum 2 ^^ at temperatures 
T < 1 K and also in recent experiments on semiconduc- 
tor nanowires^ and nanotubes^ at higher temperatures 
(down to ~ 10 K). 

Experiments on the thinnest nanowires are not well 
explained by current models of independent scattering 
events at the boundary* 8 -, even though such models de- 
scribe thicker wires quite well (see figure [lj . In partic- 
ular, data on a 22 nm Si nanowire^ shows thermal con- 
ductivity scaling with temperature up to temperatures 
of order 200 K, well above the temperature where higher 
transverse modes should be occupied; linear tempera- 
ture dependence^ is expected when only the four gapless 
modes are occupied. This paper studies the effect of co- 
herence between different scattering events in quasi-one- 
dimensional systems of variable transverse dimension, us- 
ing a large-scale numerical transfer matrix approach to 
check simple analytic models. The differences between 
bulk and boundary disorder are studied ^ 10 ' 11 , and linear 
temperature dependence at high temperature is found to 
depend on having boundary disorder rather than bulk 
disorder. 

The study of quantum effects on electrons in quasi- 
one-dimensional systems has developed steadily since the 
observation in 1988 of quantized electrical conductance 
in ballistic quantum point contacts^. The electrical con- 
ductance as a function of gate voltage shows plateaus at 
multiples of G e = 2e 2 /h, where h is Planck's constant 
and e the electron charge, and the factor of 2 results 



from spin. Hence G e is a universal conductance quan- 
tum for fermions and insensitive to material properties. 
The corresponding thermal conductance quantum has re- 
cently been measured in low-temperature heat transport 
through small insulating structures^: at the lowest en- 
ergies the only modes accessible in this system are four 
phonon modes, and the observed thermal conductance at 
the lowest temperatures is 



G = Ag = 4 
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where T is temperature and ks the Boltzmann constant. 

This paper studies how scattering and localization by 
disorder appear in the thermal conductance of long quasi- 
one-dimensional structures. Previous work on coher- 
ent scattering in phonon transport has concentated on 
the strictly one-dimensional limit 13 ; as explained below, 
there are fundamental aspects of phonon confinement 
even in the clean case that are not well captured by 
a purely ID chain. We find several important differer- 
ences between scattering and localization effects in ID 
phonon transport and the corresponding effects in elec- 
trical transport. As in the electronic problem, transport 
in quasi-lD systems can be ballistic, diffusive, or local- 
ized, but unlike the electronic problem, the lowest four 
harmonic modes are protected from scattering as T — ► 
by fundamental properties of rigid bodies (in the same 
way as the three acoustic modes of bulk solids). In order 
to access very large system sizes, many of our numer- 
ical calculations are carried out on in two rather than 
three spatial dimensions; however, two-dimensional sys- 
tems already include some important features absent in 
one dimension, e.g., different dispersion relations for the 
lowest-lying modes (u> ~ k for torsional and longitudi- 
nal, u! ~ k 2 for transverse). This paper presents two 
separate approaches to understanding the thermal con- 
ductance of these thin nanorods. Firstly, we use a scalar 
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model to study the mean free path of different phonon 
modes. We find that the gapless mode (which has no 
transverse wavevector) has a mean free path significantly 
longer than the length of the system; when combined 
with the contribution of the other modes, the resulting 
thermal conductance for very thin rods (of order 20 run) 
is approximately linear. 

We also obtain the thermal conductance at low temper- 
ature by numerical transfer-matrix calculations on sys- 
tems of small transverse dimension and an analytic the- 
ory that makes standard mesoscopic assumptions about 
the nature of scattering. Calculations are done using the 
Landauer formalism valid at low temperature o 2 ' 14 i 15 ' 16 . 
The analytic results for boundary scattering build on the 
scalar approximation introduced in Ref. 1 1 71 - The coher- 
ent multiple scattering discussed in this paper is much 
more important in long nanowires and nanotubes, where 
the ratio of length to width may be 100:1, than in the 
suspended membrane devices studied in Ref. Q], where 
nonuniform width and individual scattering effects have 
been shown to explain the experimental data&i£. 

The results obtained here suggest that the experimen- 
tally observed dip in G/T at moderate temperatures, dis- 
cussed in detail for the experimental geometry of Ref. 1 
in Ref. HI, will appear generically in one-dimensional sys- 
tems as a consequence of derealization of phonons at 
low energy and 1D-3D dimensional crossover at high en- 
ergy. Our theoretical predictions are also compared to ex- 
periments on semiconductor nanowires at relatively high 
temperature as a function of diameter—. We find that, 
while the linear temperature dependence on the small- 
est nanowire in this experiment is a robust phenomenon, 
the coefficient does not need to be close to 4go, although 
it happens to be remarkably close to this value in that 
experiment. 

We use the term "phonon" to describe both long- 
wavelength bulk modes of the nanowire and conventional 
short- wavelength phonons as there is no sharp distinction 
between these two limits. Note that phonon localiza- 
tion by strong local scattering can occur also in d > 1 
(a review is Ref. [Lsl ). but considerably more disorder is 
required. The following section discusses the basic the- 
ory of mesoscopic thermal transport in harmonic quasi- 
one-dimensional systems: localization is described by the 
DMPK model with appropriate modifications. In section 
III we argue for the relevance of Anderson localization 
and mode inequivalence for thin nanowires. Section IV 
contains numerical studies of a 2D model to understand 
the key differences between phonon and electron meso- 
scopics and test the DMPK model, and Section V ap- 
plies this model to three-dimensional nanowires at low 
temperature. 
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FIG. 1: The thermal conductance as a function of tem- 
perature for various wire diameters. The fit (solid lines) 
matches the data to a Boltzmann equation treatment includ- 
ing impurity, surface and umklapp scattering. Note that the 
22nm wire is anomalous. This data is taken from Li et al., 
Appl. Phys. Lett, 83 2934 (2003) (Ref.@); the fit (except for 
d = 22) is from N. Mingo, Phys. Rev. B, 68 113308 (2003) 
(Ref. 3); figure courtesy A. Majumdar. 



thermal current in a wire from a reservoir at temperature 
T + 6T to one at temperature T is the difference between 
that for the phonons travelling to the right minus that 
for the phonons travelling to the left: 

f dk 1 

J = E / ^^( k >( k \n^k)/k B{ T +S T) _ 



E 



dk 1 

— Tkvi(k)v(k) elUi>i(k)/kgT _ - Ti(u). (2) 



Here the summation index i runs over all propagating 
phonon modes and Ti(uj) is a transmission probability of 
phonons in mode i at frequency w, defined more precisely 
in terms of the heat flux below. The boundary condi- 
tion is that a possibly disordered segment of finite length 
is connected at both ends to infinite reservoirs with no 
disorder. Note that the distribution in the center of the 
wire is not a thermal distribution, and has no well-defined 
temperature: the Landauer approach hence does not ap- 
ply if dissipation within the wire leads to thermaliza- 
tio n 19 i 20 . The wire is assumed to be perfectly harmonic 
unless otherwise stated; reduction of thermal transport 
by anharmonic scattering is exponentially damped below 
the Debye temperature.— 

The thermal conductance G of a nanowire with per- 
fect transmission (7^ = 1) of harmonic vibrational modes 
from a thermal distribution at temperature T is G = 
J/ST, or 



II. MESOSCOPIC THERMAL TRANSPORT 



^huj/k B T 
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duj. (3) 



Our approach is based on the thermal analogue of the 
Landauer formalism for electronic transport. The net 



Here N(lo) is the number of propagating modes at fre- 
quency lu and ks is Boltzmann's constant. For a 3D 
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elastic rod, such as a semiconductor nanowire, there are 
four modes that survive down to zero frequency: one tor- 
sional and one longitudinal mode, each with linear dis- 
persion u> ~ k, and two flexural models with quadratic 
dispersion uj ~ k 2 . 

Note that this conductance is finite even if there is 
perfect transmission of heat from one end of the wire to 
another: this is the thermal equivalent of the well known 
"contact resistance" in ID electronic systems. If N(ui) is 
constant, so that N(w) — N(ui — 0), we obtain in this 
ballistic limit 



G = N{0) 



7r 2 k B 2 T 
3h 



(4) 



The same formula gives the leading behavior for small 
temperature T even if N(u>) is not constant, because in 
this limit the Bosc-Einstein factor in the integral becomes 
concentrated near u> = 0. 

Elastic scattering is included in the Landauer approach 
through the "throughput" T(u>) = ^ %{oS), which gives 
the fraction of incident energy flux on the left end that 
is transmitted through to the right end (hence T(io) = 
N(u>), the number of modes, for clean systems with no 
scattering). The conductance is then 



G = — !— / T(u 
2ttH 1 



k B T 2 (e^/k B T _ x): 



■ dco. (5) 



The numbers %,{uS) for disordered ID systems become 
exponentially small once the length of the ID wire is 
larger than a "localization length" £ that is determined 
by disorder strength and wire geometry. The analytic 
model of Maynard and Akkermans 13 on a pure ID chain 
(N(w) — 1 for all allowed cu) with random masses can be 
summarized as follows: assume the throughput T(u>) is 
unity for L < £ and zero for L > £, where the localization 
length £ for the ID linear chain follows from results of 
DysonS 2 -: 
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(6) 



Here <tm 2 is the variance of the random mass distribution 
and u>d is the Debye frequency. Note that the frequency 
dependence of £ means that the same wire can be in the 
localized regime (i.e., L»() for high frequencies, while 
in the ballistic regime (i.e., L £) for low frequencies. 

There are three regimes in a quasi-one-dimensional 
(multimode) wire. Let £{oS) be the mean free path av- 
eraged over the modes that propagate with frequency u>. 
At a given oj, if N£ ^> L then 



N(u>) 



i + L/e(u) 



O{N ) 



(7) 



This equation describes ballistic behavior for £(ui) ^> L, 
and for N£ ^>> L £, Ohmic behavior in which the ther- 
mal conductance of a long wire decreases as the reciprocal 



of the length. If N£ <C L, then the phonons at frequency 
u> are effectively localized and 



T 



(8) 



where the localization length is £ = N£. The number of 
propagating modes N clearly plays an important role in 
determining the crossover from Ohmic to localized trans- 
port. 

These conclusions follow from noting that the same 
DMPK equatio n 23 ' 24 i 25 approach valid for electronic sys- 
tems is appropriate for phonons, except that the depen- 
dence of the mean free path I on frequency is modified. 
The connection between T and the observed transport 
coefficient, expressed in the Landauer formula, is also dif- 
ferent for electrons and phonons, since in the phonon case 
all frequencies with fko < k B T contribute significantly, 
while in the electronic case only frequencies within k B T 
of the Fermi level contribute significantly. 

The DMPK equation is valid as long as phonon trans- 
fer matrices on scales longer than some mean free path 
become described by the orthogonal random matrix en- 
semble for one-dimensional systems 2 ^. While we cannot 
prove the validity of random-matrix theory assumptions 
for this problem, Section [V] verifies that the Ohmic ([7|) 
and localized ((5|) limits for T correctly describe thermal 
transport in a model system with bulk impurity scatter- 
ing, while strongly localized transport in models with 
boundary scattering is not observed. Section Hvl explains 
how microscopic transmission coefficients are calculated 
via a transfer-matrix approach. In the next section we 
analyze the case of surface scattering. 



III. STRONGLY DISORDERED SURFACE 

It is typically assumed that the mean free path of a 
phonon due to surface scattering is or the order of d, the 
width of the material. In this section we show that this 
is not true for all phonon modes. Consider, for instance, 
the gapless modes (ie. those that have vanishing wave- 
vector perpendicular to the long direction of the rod.) 
It is possible to estimate the mean free path of these 
modes using the following model. Consider a thin three- 
dimensional rod of length L, of square cross-section with 
each edge of width d. Allow scalar waves to propagate 
on the strip, with Neumann boundary conditions. The 
dispersion relation is then 



to 



k 2 



(9) 



where fcj^ = \J n 1 + n z( n /d), and n x ,n y — 0,1,2,.... 

We seek to calculate the mean free path due to surface 
scattering of the n x = 0, n y = mode. 

It can be shown using the methods of Ref. [lj] that, for 
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this mode, the mean free path is given by 



4(3tt 



2\2/3 _ 



UJ D 
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2\l/3 



(h/af 



UJ D 

(10) 

where h is the mean width of the disordered surface layer, 
c is the velocity of the gapless mode, and the correlation 
length scale for the surface disorder has been set to a, 
the lattice spacing. The first term comes from scatter- 
ing back into the gapless mode, and the second into all 
other modes, and so involves a factor of the density of 
states, which scales like uj 2 . The frequency dependence 
here is what one expects for scattering from point-like 
impurities. The gapless mode, then, does not sense that 
the disorder is at the boundary: it sees the disorder as a 
collection of point-like impurities. This frequency depen- 
dence ensures that at low temperatures the mode will be 
scattered weakly, and so will "reflect" almost specularly 
from the boundary. 

This weak scattering of the gapless mode agrees with 
Rayleigh's criterion on the condition for specular reflec- 
tion from a boundary. Rayleigh considered the follow- 
ing simplified model of a disordered surface: suppose the 
clean surface is a straight line at y = 0. Let the surface 
be at either y = h/2 or y — —h/2. An incident wave 
then reflects from the disordered boundary as from two 
Bragg planes, a distance h apart. Define 9 as the angle 
the incident wave makes with the surface. The phase 
difference between different parts of the reflected wave is 
then given by 



. , „ 2hsin6 

A</» = 2?r = 2hk± 

A 



(ii) 



where k±^ is the component of the incoming wavevector 
perpendicular to the surface. If this phase difference A</> 
is of order 1 , then the reflected wave is destroyed through 
destructive interference, and the so wave must have been 
reflected diffusely. On the other hand, if A</> <C 1, the 
reflection will be specular. 

Consider now the phonon modes of the quasi-one- 
dimensional rod we considered above. One expects then 
that the modes with lower k± will scatter more weakly 
from the boundary than those with a larger value. As we 
increase k± through the value 1/h, one expects a transi- 
tion from modes with a long, frequency-dependent mean 
free path to modes with a mean free path of order d. 

We can make an ansatz for the total throughput in the 
case of strong mode inequivalence, assuming that only 
the gapless mode has the frequency dependent mean free 
path, and all others have a mean free path of d: 



i 



L/ti 



(12) 



where li — Iq for the gapless mode and equals d for all 
other modes. This formula gives the DMPK form in the 
case of mode equivalence (^ = £), and gives ballistic 
behavior if one of the mean free paths is much longer 
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FIG. 2: The throughput as a function of frequency for a sur- 
face disordered rod, as suggested by equation (|12[1 . The length 
was chosen to be 2000 nm, a typical length for a nanowire. 
Here we have assumed that only the gapless mode has a fre- 
quency dependent mean free path, and all others have a mean 
free path of order the width of the system, d. The label 
"IBM" stands for "one ballistic mode". 




T/T D 

FIG. 3: The thermal conductance associated with the 
throughput functions shown in figure [2] The scale for G(T) is 
set by the curve for one ballistic mode (IBM). We have also 
included the correction due to umklapp scattering using the 
form suggested by Mingo^, but this correction is very small 
in the temperature range of interest. 



than the length of the system. We have assumed that the 
strongly scattered modes give an Ohmic contribution. In 
detail then, the throughput for the rod is 



1 



L£^ 



d (3tt 2 ) 2 / 3 d 2 f uj 
L 



8tt 



UJ D 



(13) 



This is plotted for various widths in figure [2] This func- 
tion can also be numerically integrated using equation [5] 
to find the corresponding thermal conductance as a func- 
tion of frequency, G(T). This is shown in figure [3] 

The throughput is the sum of two contributions: the 
gapless mode gives a term that scales like 1/(1 + Cuj 4 ); 
the other modes a term that scales like uj 2 . Together, 
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for a sufficiently thin rod, the dip in the throughput is 
somewhat smeared out by the exponential factors in the 
integral for G(T), with the result that the thermal con- 
ductance is approximately linear in T up to about Tq/5. 

It is also interesting to consider whether localization 
effects might be relevant. Let us assume that, at a fre- 
quency u>, there are of order N(u>) modes with a mean free 
path of d (we ignore for the moment the quasi-specular 
modes). For this given mean free path, the localization 
length becomes ~ N{to)d. For low frequencies, the 
number of modes is small, and so £ should be of order 
the width of the system, which by assumption is much 
shorter than the length L. Hence all the low frequency 
modes should be localized. As we increase the frequency, 
the number of modes increases like 



N(w) 



vjd 



(14) 



and so at some value we will have N(uj*)ad — L, or 



•jj- t . 



' L/a 



0.24wc 



/ 22nm\ 



3/2 



V d 



(15) 



where we've set a = 0.543 nm, the lattice constant of 
silicon, L — 2 /im, a typical value for the length of a 
nanowire, and 22 nm was chosen as a comparison scale, 
since this is the width of the wire with anomalous thermal 
conductivity as observed by Li et al^. The frequency scale 
cj* corresponds to a temperature scale of 



T, ~ 0.247b 



22nm 



3/2 



(16) 



Under the assumption then these phonons have a mean 
free path due to boundary scattering of order d, one 
should observe an exponentially small contribution to 
thermal transport at temperatures less than T*. 

These considerations are of course greatly complicated 
by the mixing with the quasi-specular modes. The 
DMPK theory assumes that the transmitting channels 
are equivalent, and it is not known how the localization 
length changes when a quasi-ballistic mode is added to 
a number of strongly scattered modes. As an ansatz we 
approximate the quasi-ballistic mode as having decou- 
pled from the strongly scattered modes. The resulting 
throughput is 



1 



N(lo) 

1 + L/l + 1 + L/d ( 



-L/(N(u))d) 



(17) 



The throughput and corresponding thermal conduc- 
tances are shown in figures H] and Localization sig- 
nificantly reduces the thermal conductance below the in- 
coherent limit. 

A full understanding would require a theory that in- 
corporates inequivalent channels, something that is not 
available at this time. Our numerical work described 
later sheds much light on this problem. 
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20nm (with LOG 



T/T D 

FIG. 4: The throughput as a function of frequency for a 
surface-disordered rod, to illustrate the relevance of localiza- 
tion, based on the formula 1 171 in the text. As before, "IBM" 
stands for "one ballistic mode" ; the "with loc" curve is based 
on the formula [17] in the text, and the "without loc" curve is 
identical to that of figure [21 




T/T D 

FIG. 5: The thermal conductance associated with the 
throughput functions shown in figure 2] Again, the scale for 
G(T) here is set by the one ballistic mode (IBM). 



Note that in the high-frequency limit where it is pos- 
sible to make a wave packet of phonons much smaller 
than the width of the wire, and this wave packet scatters 
like a particle, boundary scattering allows wave packets 
directed along the rod axis to pass through: this im- 
plies that true exponential length dependence is not ex- 
pected with boundary scattering, even when most ther- 
mal phonon modes are localized. 



IV. DYNAMICAL TRANSFER MATRIX 
APPROACH 

First recall the definition of scattering and transfer ma- 
trices for waves. Let the amplitude in mode i that is in- 
coming from the left be af, incoming from the right be 
df, outgoing to the left be and outgoing to the right 
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be Cj. The index i ranges from 1 to N(u). 



3 i(wt — fcja;) 



(18) 



Then the S-matrix maps the incoming amplitudes to the 
outgoing amplitudes 

/ : \ / : \ 



(19) 



di 

\ I \ ■ ) 

The matrix elements are given by 



r a 

t r' 



(20) 



where r, t' , t and r' are all N(u) x N(v) matrices. The 
matrix elements have the following interpretation: 
is the amplitude for a left moving wave in mode i to 
be reflected to mode j; r\j is the amplitude for a right 
moving wave in mode i to be reflected into mode j; Uj 
is the amplitude for a left moving wave in mode i to be 
into mode j; and t\j is the amplitude for a right moving 
wave in mode i to be transmitted into mode j. 

For our purposes, the transfer matrix R will prove more 
useful. It is defined by 

/ : \ / : \ 



R 



(21) 



/ 



It has the property that if there are two regions in series 
labelled and 1, then the total transfer matrix is given 
by R — RiRq. It is not difficult to express R in terms of 
the transmission and reflection amplitudes: 



R = 



t 



r't'- 1 
-t'- x r 



r r't'- 1 
t'- 1 



(22) 



The total transmission coefficient is obtained by tak- 
ing the 22 element of R\Ro, and inverting it: t' 0+1 = 
t' (l— rirp) -1 ^. This formula has a simple interpretation: 

write it as t' 0+1 = io*'i+to r i r o*i + *o r *i r o r i r d*i H • ^ is 

evident that the total amplitude for transmission (from 
right to left) is given by the amplitude to get through 
the second obstacle and then the first obstacle; plus the 
amplitude to get through the second obstacle, bounce off 



the first obstacle, bounce off the second obstacle, and 
then get through the first, ad infinitum. 

Thermal transport by vibrational modes in a multi- 
mode wire requires a few modifications to the above 
equations: starting from coupled linear equations of mo- 
tion for many ionic coordinates, we derive the form of 
the dynamical transfer matrix R that describes a longi- 
tudinal step along the wire. As an example, first con- 
sider a purely one-dimensional chain of balls connected 
by springs of constant strength k and mass m. Let Uj 
be the deviation, assumed small, of the position of ball 
i from its equilibrium location. The equation of motion 
for site i in the chain is 



mui = k(ui + i — Ui) + k(u 



i-i 



Ui). 



(23) 



The equation of motion at frequency ui can be written as 
follows: 



Ui+i \ _ / 2 — mui 2 /k — 1 \ / Ui 
m J ~ \ 1 J \Ui-i 



(24) 



Hence the matrix that corresponds to adding a single site 
in this coordinate basis is 



M = 



2 — m,L0 2 /k 



1 







(25) 



Now we convert this matrix to the basis of cigenmodes 
in order to obtain the transfer matrix R defined above. 

An eigenmode of wavevector k x satisfies (here a is the 
lattice spacing) 



M 



Ui 



A 



Ui 
Ui-l 



X = e 



I- 1\, j- (.1 



(26) 



M is a real matrix and its eigenvalues form a time- 
reversed pair: they are either complex conjugates of each 
other (propagating modes in opposite directions) or re- 
ciprocal to each other (evanescent modes in opposite di- 
rections). The eigenvalues are 



A 



2k 



± 



4k /m 



2k 



(27) 



describing propagating modes with k x real if mui 2 < 4fc 
and evanescent modes with k x imaginary if raw 2 > 4. 
Consider the regime of propagating modes: in this basis 
the transfer matrix becomes just 



R 



A 
A* 



which gives a unitary S-matrix since 

S' - I A 

A* 



(28) 



(29) 



and |A| = 1. 

The same procedure applies in a muitmiode wire ex- 
cept that now r,t,r' ,t' are matrices. There is one techni- 
cal difference that appears between the vibrational case 
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and the standard electronic case. In the vibrational case, 
it is the energy flux rather than the probability flux that 
satisfies a continuity equation. Hence the basis in which 
the S matrix is (almost) unitary is one in which incom- 
ing and outgoing excitations are normalized to unit en- 
ergy flux, rather than unit probability flux. Actually the 
whole S matrix is not unitary, but only the portion made 
up of modes with nonzero energy flu?*2£. 

Note that the transfer matrix approach used in this 
paper to make contact with mesoscopic electron physics 
works at a fixed frequency u) (equivalent to fixed energy in 
the electronic case) . The extension of the above to an ex- 
tended ID wire with random masses or spring constants 
is as follows: first the eigenmodes in the absence of dis- 
order are obtained from the transfer matrix above. Then 
the transmission coefficients for these modes through a 
disordered system are obtained by combining many indi- 
vidual segments (with different transfer matrices). After 
the transmission coefficients for a fixed disorder realiza- 
tion are obtained as a function of to, a final integration 
step gives the thermal conductance for that realization 
as a function of temperature. (Note that the first step of 
this approach, obtaining the eigenmode spectrum with- 
out disorder, is different computationally from typical 
dynamical matrix calculations in which k x is fixed and 
the frequencies lu of different phonon bands are obtained 
but obtains the same spectrum of propagating modes.) 

In the multimode case, the transfer matrix R describes 
how incident and outgoing modes at the left edge of a 
region are matched to incident and outgoing modes at 
the right edge of the region. Once there is disorder in 
the system, energy flux incident in mode i will be dis- 
tributed over reflected modes, other transmitted modes, 
and the original mode, with conservation of energy flux. 
The Landauer expression for the thermal conductance 
in Section II contains a sum over %; % is defined as 
the fraction of the incident energy flux in mode i that 
reaches the other end of the sample. In the next section 
we show numerical results on the thermal conductance of 
two-dimensional disordered strips. 



V. THERMAL CONDUCTANCE OF 
QUASI-ONE-DIMENSIONAL SYSTEMS: 
NUMERICAL RESULTS 

A. Two-dimensional strip 

This section applies the transfer matrix approach out- 
lined in the previous section to a simple harmonic "ball- 
and-spring" model system with ions moving in two di- 
mensions. Considerable previous work^l has shown in 
the absence of disorder that similar models, possibly 
including next-neighbor and bond-angle terms, can de- 
scribe the phonon spectra of quasi-one-dimensional sys- 
tems (including carbon nanotubes) with ions moving in 
three dimensions. The justification for studying a model 
system with only two-dimensional motion is that already 




FIG. 6: The acoustic longitudinal mode of a strip of triangular 
lattice. 




FIG. 7: The acoustic flexural mode of a strip of triangular 
lattice. 

in two dimensions both linearly and quadratically dis- 
persing modes are present, as in 3D, and the analytic 
predictions in section II can be checked comprehensively 
by going to very large wire sizes. The following section 
fSection lVBp sets out experimental predictions for semi- 
conductor nanowires in 3D based on the analytic picture 
outlined in Section II and checked below. 

The existence of gapless modes in a vibrational system 
is related to the Euclidean motions of a rigid body: the 
gapless vibrational modes are those that become symme- 
tries, typically translations, in the long- wavelength limit. 
There are four gapless modes in a rod, rather than three 
as in a bulk solid, because rotation around the rod axis 
guarantees a gapless "torsional" mode in addition to flex- 
ural and compressional modes that become translations 
along the three axes. For a rigid wire moving in two di- 
mensions, there are two gapless modes: using x to denote 
the rod axis, there is a linearly dispersing longitudinal 
mode with displacements along £ in the long-wavelength 
limit, and a quadratically dispersing flexural mode with 
displacements along y. Additional modes begin to prop- 
agate as the frequency is increased. 

The simplest rigid body in two dimensions is the trian- 
gular lattice with springs connecting lattice points. (The 
square lattice has vanishing rigidity in 2D if the potentials 
are purely length-dependent nearest-neighbor springs.) 
The longitudinal mode for this system is shown in figure 
[5] and the flexural mode is shown in figure [7] 

A simple model for the number of propagating modes 
agrees qualitatively with the numerical results for this 
quantity. The system has two modes as uj — > — the Ion- 
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FIG. 8: The number of modes N(u) for a strip of width 32 
sites for u < ujd- N(u>) is approximately linear in w for this 
2D system, justifying the form for N(uj) used in the text. 



gitudinal and flexural acoustic modes. As we increase uj 
to ~ Trv/2d, where d is the width, we expect to encounter 
modes that have a finite wavelength in the transverse di- 
rection. If we assume, in a simplified picture, that at 
each uj = nwv/2d we add 2 more modes to N(uj), and 
that every mode has the same velocity, then the number 
of modes is approximately given by 



N(uj) 



4d 



(30) 



For uj > lod, where ojd is the Debye frequency, the low- 
est lying modes are no longer propagating, and so one 
expects the number of modes to decrease, with approxi- 
mately the same slope. 

In what follows we will assume that N(lj) has the form 



,w s „ du 
N(oj) = 2 + c 

0J D 



(31) 



where c is a constant of order unity, and d, the width, is 
now measured in units of the lattice spacing. 



1. Bulk disordered 

Consider the case where the isotopic disorder extends 
throughout the material (in a later section we will con- 
strain the disorder to be at the edge of the strip). The 
rate at which waves scatter from point-like defects is de- 
termined by the Rayleigh-Klemens formula 28 : 



,2, ,2 



,V-1 



V- 



(32) 



where T> refers here to the dimension of the space, and 
<7 2 is the fluctuation in a single mass. The DMPK for- 
malism suggests a formula for the total throughput of a 
quasi-one-dimensional system in the ballistic and Ohmic 
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FIG. 9: The throughput as a function of frequency for a two- 
dimensional strip of width 32 sites. The various lengths are 
shown in the legend. The throughput shown has been aver- 
aged over 40 implementations of the disorder, as described in 
the text. 



regimes: 



T(uj) 



N(oj) 



1 + L/£[uj) 



(33) 



As the length of the system L extends belong the mean 
free path £, the throughput (and hence the contribution 
of that frequency to the thermal conductance) scales like 
1/L; this is a feature characteristic of Ohmic systems. 

From the Rayleigh-Klemens formula, one anticipates a 
mean free path of 



1 



1 



i(uj) 



a 2 uj 2 N(uj) 



(34) 



In the Ohmic limit this suggests a length and frequency 
scaling in the throughput of 



T(uj) 



N(uj)£(uj) N(u) 



1 



1 



L 



L uj 2 N{uj) Luj 2 



(35) 



As an example of the throughput of a bulk disordered 
system, see figures [9l [TUl The strip chosen was of width 
32 sites, and of varying lengths as shown in the legend. In 
each case, the probability that a given site was replaced 
by an disordered site is 1/32. Here the disordered sites 
were chosen to have mass 1 + a/2, whereas the original 
masses are 1. 

The general shape of the curves in figure^is consistent 
with the DMPK formula for the throughput, with the fre- 
quency dependent mean free path given by the Rayleigh- 
Klemens formula. For small frequencies, the phonon scat- 
tering rate is small, and so the total throughput tracks 
the number of modes. As the frequency increases, the 
phonons are more strongly scattered; the mean free path 
decreases faster than the density of states increases, and 
so the throughput decreases as with increasing frequency. 
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FIG. 10: This figure shows the same data as figure [9] here 
though the throughput has been multiplied by the length ol 
the system, to make the 1/L dependence more clear. The 
curve shown is a fit to A/lj 2 ' 6 . 
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FIG. 11: Throughput as a function of frequency for a bulk 
disordered strip of width 32 sites and various lengths. The 
alloy fraction here is 0.1. As the throughput passes through 
1, it tends to zero faster than any power. 



The dependence of the throughput on L and uj is con- 
sidered in figure [TOl In this figure the throughput multi- 
plied by the length is plotted as a function of frequency. 
The throughput clearly scales like 1/L in the frequency 
regime shown. The line shown in the figure is of the form 
A/u> a , where A and a are fit to the data. 

The throughput for the localized regime is illustrated 
in figure mi As the Ohmic throughput, T ~ Ni/L de- 
creases to a value less than 1, the throughput decreases 
faster than any power law. 



2. Edge disordered 

We consider in this section the throughput of the two 
dimensional strip with only the boundary of the strip 
disordered. Recall that in the case of bulk disorder, the 
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FIG. 12: Throughput as a function of frequency. The strip is 
edge- disordered and of width 32 sites. 
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FIG. 13: Throughput of a strip of width 32 sites multiplied 
by L ' 55 as a function of frequency for the case where the 
disordered is only at the edge. 



throughput shows a strong frequency dependence — es- 
sentially because the mean free path is set by Rayleigh- 
Klemens formula for scattering from point-like defects. 
In contrast, when the disorder is present only at the edge, 
a more appropriate model is that the phonons are scat- 
tered, specularly or diffusively, at the boundary. 

The throughput is shown as a function of frequency for 
various lengths in figures [12j HSl and ITU The width of the 
strip in figures [T21 and IT51 is 32 sites, and in figure [TH is 48. 
The throughput was calculated by disordering a region 
of thickness 1 site at both edges of the strip. The mass of 
these edge sites was chosen at random to be either 1 (the 
mass in the clean case) or 1 + \/2, each with probability 
1/2. The throughput was then calculated numerically 
as described earlier. The figures show the throughput 
averaged over 40 implementations of the disorder. 

From the figures one can see that once the phonon fre- 
quency is no longer in the ballistic regime, the through- 
put levels off to approximately a constant. A throughput 
that is constant as a function of frequency can be seen 
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FIG. 14: Throughput of a strip of width 48 sites multiplied 
by L ' 45 as a function of frequency for the case where the 
disordered is only at the edge. 




FIG. 15: Mode-resolved throughput of a strip of width 32 
sites as a function of frequency for the case where the disor- 
dered is distributed throughout the bulk of the material. The 
length of the strip is 1600 sites. 



experimentally as a thermal conductance that is linear 
in temperature: if the throughput is equal to say c, then 
the thermal conductance is ego, where go is the quan- 
tum of thermal conductance. The throughput has been 
rescaled to remove the length dependence. In both cases 
the throughput scales like 1/L a , where a is less than 1: 
the system is not in the Ohmic regime. 

To illustrate the reason for the frequency independent 
throughput, we show the mode- resolved throughput in 
figure [TBI In the frequency range in which it is constant, 
the throughput is dominated by two modes, which hap- 
pen to have a small displacement at the boundary. Their 
contribution is approximately constant. 

When the disorder is exclusively at the boundary, the 
boundary displacement of the different modes will ob- 
viously determine how strongly they are scattered — a 
mode with vanishing boundary displacement will pass 
through the material ballistically. 
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FIG. 16: Mode-resolved throughput of a strip of width 32 
sites as a function of frequency for the case where the dis- 
ordered is distributed only at the edge of the material. The 
length of the strip is 1600 sites. 



In addition, these boundary displacements can vary 
dramatically for different modes at the same frequency. 
For instance, when u) is such that lu > v/d, the un- 
gapped longitudinal and flexural modes become surface 
waves. This occurs because a surface wave must decay 
exponentially into the bulk. Its dispersion must then 
be w ~ v \J^i — k^.- Therefore for a given ui, these sur- 
face modes have the largest k z among all the propagating 
modes: this identifies them as a linear combination of the 
flexural and longitudinal modes. 

By contrast, the first set of gapped modes tend to 
have an anomalously small displacement at the boundary. 
This is illustrated in table|TJ For a frequency of ui = 0.57, 
and a width of 32 sites, the k z values of the various prop- 
agating modes along with the amount of displacement at 
the boundary is shown. By boundary displacement we 
mean the displacement summed over the two sites at the 
edge of the strip. The modes here are normalized to 1/2: 
a mode with boundary displacement squared 1/2 would 
exist exclusively at the top and bottom sites of the strip. 



B. Three-dimensional rod 

In this section we consider numerical data for a three- 
dimensional rod. The system we study has 127 atoms 
in a cross-section, as shown in figure [TT] The atoms are 
arranged in a triangular lattice, with springs connect- 
ing nearest neighbors and next-to- nearest neighbors. At 
small frequencies there are four normal modes: longitu- 
dinal, torsional and two flexural. The longitudinal and 
torsional modes disperse linearly, while the flexural dis- 
perse quadratically. 

The rods are disordered in the same way as described 
for the strip. In the clean case each atom has mass 1; for 
a disordered system each atom is replaced by an atom of 
mass 1 + V2 with a certain probability, the alloy fraction. 
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TABLE I: Boundary displacements for the 14 propagating 
modes at frequency ui = 0.57 for a strip of width 32 sites. 
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FIG. 18: The number of propagating modes for the rod stud- 
ied as a function of frequency. 
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FIG. 19: Disorder-averaged throughput for a rod of length 
of length 100 sites as a function of frequency. Both the bulk 
and boundary disordered cases are shown. 



FIG. 17: A cross-section of the nanowires considered. There 
are in fact springs connecting next-to-nearest neighbors as 
well as nearest neighbors, although they are not shown here. 

As before, we will consider both bulk (i.e., each atom in 
a given cross-section can be replaced) and boundary (i.e., 
only those atoms at the outer edge of the cross-section 
can be replaced) disorder. 

We find the throughput using the transfer matrix 
method. In figure 1191 we show the disorder-averaged 
throughput for a rod of length 100 sites. Both bulk and 
boundary disordered cases are shown. The alloy fraction 
for the bulk case was chosen to be 0.1; for the boundary 
case it was chosen to be 12.7/36: this ensures that the 
number of disordered atoms in a cross-section is the same 
for the bulk and boundary disordered cases. 

The data shows that the throughput in the case of 
boundary scattered rods is noticeably less than the bulk- 
disordered case. This is because phonon modes in a rod 



tend to have more displacement at the edge than at the 
center, since at the edge the masses have greater freedom 
to vibrated 

One feature which is not apparent is the presence of 
modes which have a displacement at the boundary orders 
of magnitude smaller than that of the other modes (recall 
that these modes dominated the throughput of the edge- 
disordered strip) . This may be due to the relatively small 
diameter of the rod considered here: in the data shown, 
we have not in fact reached the limit where the surface 
waves have formed. 

Finally, it is interesting to note that the reflection and 
transmission matrices for a short system can be used to 
find the incoherent limiting conductance. If the trans- 
mission and reflection matrices for a system of length SL 
are T and i?, then the transmission and reflection ma- 
trices for a system of length nSL (denoted T n and R n ) 
are 

T n =f{l-R n ^ 1 R)- 1 T n _ l (36) 
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Rn = Rn-1 + T n _xR(l — Rn-iR) 1 T„_i (37) 

This is of course an approximate description of the trans- 
port of waves: it is analagous to adding the probabilities 
rather than amplitudes, and so it neglects interference 
effects. However, a physical system will have a phase- 
breaking length scale; for systems longer than this, the 
transport will be incoherent. Physically one would ex- 
pect this phase-breaking scale to depend on frequency 
and temperature. This may prove useful to isolating the 
incoherent signal in future studies. 

VI. CONCLUSION 

In conclusion, we have studied the thermal conduc- 
tance of insulating nanowires at low temperatures. For 
the thinnest wires, with diameters of order 20 nm, it 
was argued that frequency-dependent surface scattering 
is relevant. In particular, modes with a small k± (ie. the 
component of the wavevector in the direction perpen- 
dicular to the long direction of the wire) should scatter 
almost specularly from the boundary, whereas those with 
a larger k±_ presumably have a mean free path of order 
d, the width of the wire. The resulting thermal conduc- 
tance is shown in figure [3J for a number of diameters. 
For a width d — 20 nm, the thermal conductance is al- 
most equal to that of a single ballistic mode. This can 
be understood as follows. There are two contributions 
to the thermal conductance: that of the quasi-specular 
modes, whose number is fixed but whose contribution 
decreases with increasing frequency; and the diffusively 
scattered modes, whose mean free path is fixed but whose 
number increases with frequency. These decreasing and 
increasing terms together give a approximately constant 
throughput for low frequencies (figure [2]). 

These considerations also suggest strongly that local- 
ization effects are relevant to transport in nanorods of 
transverse dimension d < 20 nm. This was demonstrated 
by applying the DMPK formulism, which at certain fre- 
quencies is not valid — DMPK requires mode equiva- 
lence, and at small frequencies the quasi-specular mode 
has a mean free path much longer than the other modes. 
Nevertheless, it would be of great interest to explore ex- 
perimentally localization phenomena in these nanowires. 
For instance, a different temperature dependence in the 



thermal conductance and heat capacity, with the heat 
capacity measured through a bulk contact to a reservoir, 
would be strong evidence that some of the phonon modes 
are localized. 

In order to understand the relevance of coherence ef- 
fects, we have also simulated phonon transport using 
a dynamical transfer matrix method. Our simulations 
show a mode inequivalence for the case of scattering from 
a disordered surface that is not present when the impuri- 
ties are placed at random throughout the strip. The to- 
tal throughput of the surface-disordered strip is approx- 
imately constant in frequency, consistent with a linear 
thermal conductance. 

If the ballistic and strongly scattered modes were 
decoupled, one would expect some modes to have a 
throughput of approximately 1, and others to have an 
exponentially small throughput. Although the mode in- 
equivalence we observe numerically is not that strong, 
the total throughput is consistent with the picture that 
some modes are quasi-ballistic, and the rest are strongly 
scattered. Also, the throughput for the surface scatter- 
ing case shows a non-Ohmic behavior approximately de- 
scribed by a power-law: T(lj) ~ l/L a with a < 1. This 
is consistent with the presence of quasi-ballistic modes. 

An additional effect brought to light by our simulations 
is the highly variable displacement that phonon modes 
can have at the boundary of the nanowire. We found 
that whereas the gapless modes acquire a large bound- 
ary displacement, other modes acquire anomalously low 
values. These were the modes that gave a quasi-ballistic 
contribution to transport. Our dynamical transfer ma- 
trix method could be applied to specific crystal structures 
and materials of interest by using parameters from full 
ab initio simulations of nanowires, including boundary 
effects. 
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